SUBROUTINE CALC_TAU
!
! Verwendungszweck: 
!    Die zeitliche Daempfung des adaptiven Gitters kann man raemlich abhaengig von der Loesung gestalten.
!
!    Dabei stellt sich die prinzipielle Frage ob man zur Berechnung von 'tau' einmal zu beginn der
!    ausrechnet (beim updaten von XA) und dann waehrend der newton-iteration gleichlaesst
!     -- oder ob man 'tau' nach jedem iterationsschritt updated.
!
!    Mir ist nicht klar welches die bessere Methode ist:  tau und der entsprechende aktuelle Variablensatz X
!    bilden nicht (so wie z.b. X unt f_edd) ein zusammenpassendes paar. 
!    Die zeitliche daempfung haengt von tau/tst ab - und tst ist ja auch kein 'bestandteil' des variablensatzes.
!
!
!    Das hier ist nur ein Dummy-file, mal sehen ab man das fuer Scheiben brauchen kann !xxx debugging
!
!

      use gridvar, only  : tau, taurel


      implicit none


!-----------------------------------------------------------------------
!    Berechnung des zeitlichen Gitterdaempfungsparameters
!-----------------------------------------------------------------------

      tau(:) = taurel


END SUBROUTINE CALC_TAU
